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GRAVITATIONAL FIELD MODELS FOR THE EARTH 

(GEM 1 & 2) 

ABSTRACT 


Two models of the earth’s gravitational field have been com- 
puted at Goddard Space Flight Center. The first, Goddard Earth 
Model 1 (GEM 1), has been derived from satellite tracking data. 
The second, Goddard Earth Model 2 (GEM 2), has been derived 
from a combination of satellite tracking and surface gravimetric 
data. The geopotential models are represented in spherical 
harmonics complete to degree and order 16 for the combined solu- 
tion and complete to degree and order 12 for the satellite solution. 
Both solutions include zonal terms to degree 21 and related sat- 
ellite resonant coefficients to degree 22. The satellite data con- 
sisted primarily of optical data processed on 300 weekly orbital 
arcs for 25 close earth satellites. Surface gravity data were 
employed in the form of 5° x 5° mean free-air gravity anomalies 
providing about 70% world coverage. Station locations were 
obtained for 46 tracking sites by combining electronic, laser, and 
additional optical tracking data with the above satellite data. Anal- 
ysis of the radial positions of these stations and a value of mean 
gravity on the geoid indicated a mean equatorial radius for the 
Earth of about 6378145 meters. Results of geopotential tests on 
satellite data not used in the solution show that better agreement 
was obtained with the GEM 1 and GEM 2 models than with the 1969 
Smithsonian Standard Earth II model. 
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GRAVITATIONAL FIELD MODELS FOR THE EARTH 

(GEM 1 & 2) 


I. INTRODUCTION 

The establishment of an accurate Earth model (geometric and gravimetric) 
is an essential requirement of any earth physics program (1). Progress in this 
are:, has been made at the Smithsonian Astrophysical Observatory using optical 
tracking data (mainly) and analytic techniques in the computation of orbits' 2 ’ 3 ). 
Work at Goddard Space Flight Center has proceeded with the object of extending 
the use of this optical data and also to the inclusion of highly accurate and dense 
electronic satellite tracking data available since 1965 (i.e. Tranet Doppler, 
Goddard Range and Range-Rate, and laser). Numerical integration, providing 
a precise solution for modeled forces, has been used extensively at GSFC for 
the computation of orbits. With this approach a preliminary solution for a geo- 
potential model has already been obtained using 17 satellites and Baker-Nunn 
optical data only(^). in this solution the geopotential field was complete to 
degree and order 8 in spherical harmonics and it included station coordinates 
for 13 optical tracking sites. The combined solution reported here is com- 
plete to degree and order 16 and is based on surface gravity as well as satellite 
optical data. The station coordinate solution has also been extended to 46 
stations. In this new solution considerable improvement has been made through 
refinement in the modeling for drag and satellite resonant effects. 

Tests of the new solution are presented which demonstrate improvement 
over other solutions in satellite orbit determination and better agreement with 
ground survey data. Consistent results have been obtained for an adjusted 
scale of the reference ellipsoid from analysis of both station positions and 
surface gravimetric data. 
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n. DESCRIPTION OF SOLUTIONS 

A geopotential solution (GEM 1) has been derived from tracking data (pri- 
marily optical) on 25 close earth satellites. Spherical harmonic terms in the 
geopotential are complete to degree and order 12 with zonals extending to degree 
21 and selected satellite resonant coefficients to degree 22, Surface gravity 
data, consisting of 21,000 1° x 1° mean free-air gravity anomalies and used in 
the form of 5° x 5° mean anomalies, have been combined with the satellite data 
to derive a geopotential (GEM 2) complete to degree and order 16 with similar 
higher degree terms as in the satellite solution. 

Station locations were obtained for 46 tracking sites by combining satellite 
electronic, laser, and additional optical data with the data used in the satellite 
geopotential solution described above. Tracking systems consisted of 13 Baker- 
Nunn cameras, 23 Minitrack Optical Tracking System cameras (MOTS), 2 
lasers, 3 Goddard Range and Range-Rate (GRARR) systems, and 5 NWL Tranet 
Doppler systems. A value of mean gravity for the reference system was deter- 
mined through a simultaneous adjustment of this quantity with the geopotential 
coefficients using both the satellite and surface gravity data. 

Data Used in the Solutions 

a. Satellite Data: Data used in the satellite geopotential solution (GEM 1) 
is presented in Table 1, Approximately 300 weekly orbital arcs of optical data 
on 23 satellites formed the basis of the solution. Also included in the solution 
were 21 weekly arcs of Minitrack interferometer data on TIROS 9 and Allouette 
2 satellites, which v/ere principally used to support the determination of sat- 
ellite resonant coefficients. Additional information is presented in Table 1 on 
satellite orbit geometry and the number of orbital arcs, observations and 
observation residuals. 

In the solution for the 46 station locations, the data in the above 300 weekly 
orbital arcs (120,000 optical observations) were combined with 15 weekly GEOS-I 
and n orbital arcs of densely covered electronic and laser data (150,000 observ- 
ations) including some additional optical data. Also, some 66 one- and two-day 
arcs of GEOS-I and n flashing light optical data were combined with the above 
arcs, mainly to strengthen the solution of the MOTS stations. 

b. Surface Gravity Data: The source of most of the gravimetry data was 
the U. S. Aeronautical Chart and Information Center which provided 19,000 one- 
degree by one-degree mean free-air gravity anomalies. A further set of 2000 
mean gravity anomalies were obtained from a number of other sources. These 
data were used to form 1707 five-degree by five-degree mean gravity anomalies 
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by a straight averaging of the one-degree by one-degree mean anomalies, and 
provided a total coverage of about 70% of the earth’s surface. See Figure 1 for 
a map of surface gravity data coverage. 

Starting Values for Solutions 

The a priori values of the geopotential coefficients and station locations 
used in the solution were as follows: 

Zonals harmonics to degree 21, Kozai (1969)(^) 

Tesseral harmonics SAO S.E. I (1966)( 2 ) 

Satellite resonant coefficients GSFC preliminary analysis^) 

Station coordinates Marsh, et al. (1971) C 7 ) 

The reference ellipsoid that was adopted for the solution was the same as 
that used by SAO in the Standard Earth II (1969)( 2 ); viz 

Mean equatorial radius, a e = 6378155 meters 
Flattening, f = 1/298.255 

Product of mass of the earth and gravitational constant, 

G-M = 3.986013 x 10 14 m 3 /sec 2 . 
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III. ANALYSIS TECHNIQUES 

The motion of the satellite was obtained by numerical integration of the 
equations of motion with the following forces being modeled: the earth’s gravity 
field, atmospheric drag, solar radiation pressure, and lunar and solar gravity. 
The integration was performed in a cartesian, geocentric, inertial coordinate 
system referenced to the true equator and equinox of the epoch. The time frame 
was A1 time and modeling in the reference system included luni-solar precession 
and nutation of the earth and polar motion, the latter being obtained from the 
U. S. Naval Observatory who also provided the adjustment from UTl to Al. The 
origin of the coordinate system for the pole positions was the mean pole of 1900- 
1905. The technique of weighted least squares was applied to the satellite 
observation equations with the standard deviations presented in Table 2. 

The surface gravity data was given in terms of 1° x 1° mean gravity 
anomalies (mgal) referenced to the International System. This data was used 
as described below in terms of the reference system defined by GM, a e , f, and 
co (rotation rate of the Earth). 

Let W be the gravity potential of the Earth so that 

W = V + (1) 

where V is the gravitational potential, and <h the centrifugal potential, and 
where 


°° n / \ n 

+ S 2 W P " (sin ^ (C nm cosmX + S nm sin mX ) 
n = 0 m = 0 

1 , 

$ = -^ (r co cos < pr 

r, <p, X spherical coordinates of radial distance, latitude, and longitude, 

P™ (simp) associated Legendre polynomial of degree n, order m and argu- 
ment sin <p, C nm , S nm spherical harmonic coefficients. 

On a equipotential surface (i. e. , mean sea level), we may write W = W Q 
so that gravity g 0 , on this surface is given by 



g„ = IAW 0 | 
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The quantity g 0 is related to the international System for gravity measurement 
as follows 


go' = + ?i 

So ~ So'" 13 - 7 ( 2 ) 


where 

Yj is normal gravity on the International Ellipsoid^ 8 ) 

Agj is gravity anomaly for the International System 
-13.7 mgal is the Pottsdam correction. 

The 5° x 5° mean gravity anomalies were formed from a straight average 
of the given 1° x 1° mean anomalies for each 5° x 5° section and were repre- 
sented at the mid-point of the section in latitude and longitude. A value of g 0 
was then computed from (2) using each of the 5° x 5° mean anomalies, and the 
resultant g 0 was taken as the observed quantity. A residual, g 0 - g c , was then 
formed for the equation of condition and expressed in terms of linear adjustments 
to the starting values of the spherical harmonic coefficients,, The quantity g c 
was computed as 


_j/3W\ 2 , 1 / 3W\ 2 | ' A 

^'Kar/ rHaJl <=>> 

where the partials are obtained from (1) with the use of our starting coefficients 
and reference system parameters. A constant term Ag 0 ^ 9 ), variation of mean 
gravity from the reference system, was added to the equation of condition for 
simultaneous adjustment with the spherical harmonic coefficients. Weighted 
least squares normal equations were then formed for simultaneous adjustment 
with the satellite normal equations. 

In forming the weighted least squares normal equations the variance a 2 
of g 0 (observed gravity value corresponding to the 5° x 5° mean anomalies) 
was obtained from 


a 


2 


N + 1 


(4) 
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where 


N is the number of 1° x 1° mean anomalies within a 5° x 5° section and a^, 
the computed variance of the 21,000 1° x 1° observed anomalies, is 33 2 (mgal 2 ). 

For an unobserved 5° x 5° area, a value of zero was adopted with a variance 
of <Jq . The weighting of the gravity data was according to the formula 

cos (b 
w = 

Cr2 (5) 

where the cosine of the latitude was introduced to account for the unequal area 
covered by the 5° x 5° section. For a 5° x 5° mean anomaly derived from 25 
1° x 1° means, o - 6 mgal. 
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IV. RESULTS 

1. Solutions 

The geopotential and station coordinate solutions were described in a 
previous section, Description of Solutions. The values of the geopotential co- 
efficients, for the satellite solution (GEM 1) and the combined satellite/surface 
gravity (GEM 2) solution, are listed in Table 3. The station coordinates for the 
46 tracking stations are listed in Table 4. The ellipsoid heights for the 13 Baker- 
Nunn stations have been lowered by 15 meters to account for an error in the 
computer program for the parallactic refraction correction. The solution value 
for Ag 0 , the variation of mean gravity from the reference system, was 3.28mgal. 

2. Zonal Coefficients 

Zonal coefficients complete . j degree 21 in the GEM 2 solution are obtained 
from the combined effects of both satellite and surface gravity data. These 
zonal values are compared in Table 5 with the SAO 1966 (Ml) and 1969 Standard 
Earth n values derived by Kozai from secular and long period orbit perturbations. 
The table also contains the zonal values from our satellite only solution. Since 
weekly satellite arcs were used in our solution the satellite zonal coefficients 
are obtained from the short term zonal effects. The comparison in the table in- 
dicates remarkable agreement with SAO S.E. n in the zonal values from our sat- 
ellite only solution and even better agreement in our combined solution. 

A zonal profile of geoid height is given in Figure 2 for the GSFC satellite 
and combined solutions and for the SAO Standard Earth II solution. Very little 
difference exists except in the region of 0 to ± 20° latitude which is given on an 
enlarged scale in Figure 3 where differences of a few meters may be seen. 

Generally high correlation exists between consecutive odd and consecutive 
even zonal coefficients in satellite solutions from secular and long term orbital 
effects. Correlations, although very large (generally greater than 0.9) for 
consecutive odd or even zonals, fall off as odd or even pairs separate in degree. 
However, in our combination solution the coefficients do not exhibit these high 
correlations and have essentially become decoupled through the surface gravity 
data. 

3. Geoid Height Maps (Undulations) 

Geoid height contour maps at ten meter intervals are presented in Figure 
4 for the GSFC combined and the Standard Earth n olutions. The two solutions 
generally compare very favorably but differences o> 8 meters can be seen 
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between the relative highs and lows on the maps. An rms of 5 meters was 
estimated between the GSFC combined and the S E. II solution from a more 
detailed analysis, A maximum difference of 15 meters exists at -70° lat. and 
20° W long. , where the GSFC combination geoid is shaped differently from the 
SAO solution. 

4. Degree Variances of Gravity Anomalies 

Degree variances of gravity anomalies have been obtained, for each degree 
n = 3 to 21 , as follows 

a g 2 = t 2 (n - D 2 ± Q (Cj + S n ^)(mgal) 2 (6) 

where y is a mean value of normal gravity and C nm and S Iltn are the spherical 
harmonic coefficients in normalized form. These values are tabulated in Table 
6 for our combined and satellite only solutions and for the S.E. II solution for 
purposes of comparison. Although the total sums of the degree variances differ 
by 3 mgal 2 between the S.E. II and our combined solution, individual differences 
are as large as 8 mgal 2 . Comparisons of the solution with observed gravity a- 
nomalies are in process but are not available at the present time. 

5. Geopotential Coefficients and Standard Deviations 

Tables 3a and 3b show the geopotential coefficients obtained in the GSFC 
satellite and combination solutions. Table 7 presents an rms of geopotential 
coefficient differences with the Standard Earth II for each degree for a variety 
of solutions. Two other solutions in addition to the GSFC solutions are presented. 
These are the SAO 1966 Standard Earth I field and another SAO solution, the 
B13.1 which is similar to the Standard Earth n. It is noted that the smallest 
rms total is for the GSFC combined solution and the largest total is for the Stand- 
ard Earth I. 

The rms differences per degree between the SAO S. E. n and the GSFC sat- 
ellite solution have been used in an analysis for determining the relative weighting 
factor between the surface gravity normal equations and the satellite normal equa- 
tions . In Figure 5 the average standard deviations of the coefficients for each degree 
n are plotted for the satellite and surface gravity solutions. It might be expected 
that the two curves would cross at some point, such as at degree ten, above which the 
surface gravity would provide a relatively stronger contribution to the higher degree 
coefficients . However, the disparity between the two curves suggests the satellite 
values are unrealistic. A more realistic standard deviation of the satellite coefficients 
can probably be obtained from the rms of coefficient differences between two 
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solutions containing satellite data. This was done with the GSFC satellite solu- 
tion and the Standard Earth II and the resulting curve is shown in Figure 5. 

This curve crosses the surface gravity curve at degree 10 with a value of about 
5, 5 x 10 The satellite standard deviation at this point is 2.0 x 10' 8 so that 
the ratio xor scaling of the satellite standard deviations is 2.7. This corre- 
sponds to a relative increase of weight of 7 for the surface gravity normal 
equations. In the combined solution a weight of 5 was finally used because the 
standard deviations for the surface gravity equations, although more realistic, 
may be somewhat optimistic. 

It should be pointed out that the correlations among the non-zonal geopoten- 
tial coefficients are generally very small. Under these conditions the relative 
weight of separate normal equations in estimating a parameter for a combined 
solution is nearly inversely proportional to the ratio of the squares of the standard 
deviations associated with the separate normal equations, 

6. Comparison of Station Heights with Local Survey 

The station heights above the geoid have been obtained for 46 stations using 
our combined satellite-surface gravity geopotential and station coordinate 
solutions. These heights were then compared with the mean sea level heights 
given by local survey. The comparison is presented in Figure 6. The mean 
sea level height, MSLH, for each of the stations was obtained from the survey 
sheets in the NASA Directory of Observation Station Locations published at 
Goddard Space Flight Center. The MSLH heights are listed in Table 8 for each 
of the stations along with the corresponding geoid height, N, obtained from our 
solution. The survey values of MSLH are generally listed with an accuracy of 
a meter or better. The full station coordinates obtained in our combined 
solution are listed in Table 4 which gives the ellipsoidal height (h), latitude, 
and longitude of all 46 stations. 

These 46 stations are plotted in Figure 6 in terms of a letter signifying the 
tracking station data type, namely B - Baker-Nunn (13), M - MOTS (23), D - 
Doppler (5), G - GRARR (3), and L - Laser (2). 

The quantity plotted in Figure 6 for each of the 46 stations is 

AH = h - N - MSLH (7) 

as a function of latitude of the station. Since the geoid heights, N, and the 
station heights, h, are referenced to our ellipsoid of a e = 6378155, the zero 
value corresponds to this value of a g on tk-j plot. The average AH is close to 
-10 meters, implying a e = 6378145 meters is a better reference radius for this 
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set of stations. This value of AH = -10 meters is based upon the MOTS, laser 
GRARR, and doppler coordinates; the Baker- Nunn coordinates were excluded 
from this calculation because of the adjustment for parallactic refraction made 
to the station heights, as noted above in section 1. 


In Figure 7 a similar set of AH results were computed using the Standard 
Earth n(3) geopotential and station coordinates for 37 sites. These results 
show an average value >f about a e = 6378137. There appears to be a systematic 
difference in the sc ale implied for the Northern than for the Southern Hemisphere, 
due principally to the European stations. These stations are circled in the fig- 
ure and when excluded an average a e = 6378143 results. 

7. Adjusted Ellipsoid Scale (a e ) Based upon the Mean Value of Gravity 
Variation (Ag 0 ) 

The mean value of gravity variation Ag 0 from the reference system was 
obtained as part of a general solution as previously described. The value of 
Ag 0 = 3.28 mgal was obtained. Using the simple relation 


GM 



the variational relationship becomes, 

Agg AGM 2Aag 

IT = ~GM ~ ~ (9) 


Since the GM of the reference ellipsoid includes the atmosphere*, and the 
relative mass of the atmosphere to the earth is 0. 85 x 10' 6 and further, since 


Ag e 

8e 


Agp 

8e 


= 3.35 x 10' 


6 


where g e = 9.780291m sec' 2 , we can obtain 

a e (3.35 + 0.85) 

Aa e — - - r 

1 0 6 2 

or -13 meters. This implies that a e = 6378142 meters is a better scale for our 
reference ellipsoid. 


*GM was obtained by the Jet Propulsion Laboratory^®) from space probes and thus includes the mass of 
the atmosphere, while the free-air gravity anomalies were made at the surface of the earth and reduced to 
sea-level. 
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V. TESTS OF SOLUTIONS 


Two tests have been applied to the solutions that have been obtained so far. 
Both concern the ability of the gravity fields to represent the motion of close 
earth satellites. 

The first test was on 22 six-hour orbital arcs of the Beacon Explorer C 
(BE-C) spacecraft composed of laser range data from the Goddard experimental 
laser. The rms of fit to the data on the 22 arcs using the Standard Earth n 
gravity field varied between 0. 9 meters and 4.3 meters with a mean of about 
3 meters. With the GSFC satellite solution these rms values dropped to between 
0.6 meters and 2.0 meters with a mean of about 1.3 meters. This considerable 
improvement was obtained even though no laser data on BE-C was used in the 
solution. These rms values are shown in Figure 8. 

The second test involved the long-term behavior of INTELSAT 2-F1, Cosmos 
41 rocket and Cosmos 382 rocket. The variation of mean elements over periods 
of hundreds of days have been compared using the Standard Earth II and GSFC 
combination gravity fields. The results are shown in Table 9 from which it is 
evident that a three-fold improvement has been obtained with the GSFC com- 
bination solution. Further, it should be remembered that none of these satellites 
were used in the solutions for the gravity field. 
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VI. SUMMARY AND CONCLUSIONS 

The first comprehensive GSFC Earth Model has been obtained. This 
consists of a combined solution for the geopotential field using satellite and 
surface gravity data (GEM 2) as well as the associated satellite only solution 
(GEM 1). It also includes a set of 46 geocentric station positions. Tests and 
comparisons presented in the report yield encouraging results for these solu- 
tions, particularly in their capability for improved orbit determination as seen 
in the section on Test Results. Analysis is continuing with additional satellite 
and surface gravity data to derive a more complete geopotential model with 
spherical harmonics complete to degree 20. Satellite data from over 71 tracking 
stations are being processed for this solution. 

In the present solution good agreement has been obtained between analyses 
of surface gravity data and station height data to provide for an adjusted scale 
of the equatorial radius (a e ) of a mean earth reference ellipsoid. A value of 
a e = 6378145 is adopted. Our present solutions, data used, and results obtained 
in this report are summarized in Table 10 for convenient reference. 
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Figure 1. Surface Gravity Data 5° a 5° Mean Gravity Anomalies 






CO 

z 

o 


Z> z 
O S 

CO •" 


>- 

z 

o 


o 

CO 

H 


< X 

CO 

a: 

03 < 

_ 1x1 

Q 

LaJ Q 
Z CC 


CO 

S 

o 

o 

o 

u_ 

CO 


< 

a 

z 

< 

I— 

CO 

o 

< 


O CO 



o 

CO 


X 

o 

CM 


i 


o 


o 

r-H 

I 


(SU313IAI) 1H0I3H 


18 


GSFC COMB. 
& SAO 


-20 


i — — r 



o 

CO 

I 


o 

<T> 


O 

CO 


o 

CO 


O 


CO 

LU 

UJ " 

DC % 

O £ 

LlJ g 

Q o 

N 


O 

CO 

I 


LlJ 

a 


H* 

< 


_c 

a) 

’5 

X 

T 3 

'5 

a> 

a 

CN 

2 

3 

o> 

LL 


o 

VO 

I 


o 

<J> 

I 







< -30 


-60 


-90 

-180 -120 -60 0 60 


LONGITUDE DEGREES 

STANDARD EARTH S. SOLUTION 









DEGREE n 

Figure 5. Average Standard Deviation of Coefficients per Degree n of Satellite and Gravimetiy Solutions 


21 


<SU3i3W) H1SW N-M * HV 


o 

co 



c 

o 


C/5 


CC 

O 


LU 


o 


H 

O 

£ 

< 

CO 


LU 

O 

D 


< 


_3 

O 

to 

o 

LL 

to 

(J 


> 

0) 

> 


^3 

CO 


~o 

'o 

Q> 

a 


<u 

> 

O 

_Q 

< 


-C 

o> 


a> 

X 


c 

o 


co 


CO 


CD 

£ 

Oi 

ii. 


(SH313W) H tSW*N M=H7 


22 



AH = h-N-MSLH (METERS) 


o 

CSl 


o 

m 


o 



23 


LATITUDE OF STATION (DEGREES) 

Figure 7. Station Height Above Geoid vs Survey (SAO Solution) 





Using BE 





O v? 

tu 

tt’ 

m O' 

•0 

>o 

e 

m 

*- m 


fs 

<r C 

f.. 

K N 

w 

V 

o 

n 

c» 

O' 

N -• 

N 

t*l 

O' e> 

f\J 

< 

m 

IT 

<.* N 

1*} 

•J 

'O O 


a <t 

*• 

* 

OJ 

** 

(M 

<V 

V* «“» 


« 

*4 O 

w 

- 


Ci 

-- 

w 

vv 

-* »-l 

- 

«* *-> 


im 


104 

IM 

M i 


CQ 

O a « « n c <t ij « a 4 4 'f. fi 

n ftjoGNwowo«woinij , ^XK(?-*’ 0 ©^^' ft W '0 2 

^ »* -« D < 1 * s ft; N Ot O ^ *V •* « N M ” 


* * # * 

<j a‘ sfl o c o o o k ci *11*. o i/i is o o » « e sr c ? « o 

~» ~ in n<o « cm*-. n -* -* c 

<* w ~ 


■» 4 > SC 1 N 4 04 snOOHP^«SCNO*«'fl»| 
n — <M 4 “ r\i oj ~ ** I 


nC'oam-oh^c'iftNnfv^ r><fi«&‘cnON' 0 * 

-O 0 'N«C'WOrthf»®'f<c-ICNN 8 »O* 4 O 

Ci pj « « ^ c: p) n o n tv n n « <i x m f) w n i ^ m « 


n Mi) « 6 OOCMCC)NS!flinnO 3 )SO 0 lCNft 0 < 

-«SO<f"J-C-C«MS 6 l - — C -* <» W O <J CWN 
in oe-t-*Oox<j®N**^«c ufli{>oft:!M'JO(0 , a-* 
0 t^P*w<flAU 5 ’J)C‘alO- 4 'Cv'^® 9 ^ D OIIMC tf) (P 


m 

« N'fl xi/ltt»*NNc 05 ftiJ^^<?f)JiflOw ,;, 1 }tO o 

cv.dejomofi^u i ^K 6 ^inNOK 4 nN«>iM^ 4 to « 

♦ ♦*p«s^)N(j«)i!)dnji:o-«o*»Noe#«'e'C» 8 

CJ^OO^^OOOOOOOOOOO ooooo****.— o 

OOOCOOOOOOOOOO^OOOOOOOOOO $ 


^ H «10 5 >{g / 


c s at 

Mp) I u 

i ; ~ o' *- 

! y U if « < I ^ 


in <n 
« m n 

I I I 
U j Q 

« s * a « x 

I I I < < < 

01 ? < '/» ? ? 3 

1 3 2 <f i> J ; 3 

^ l n £ 2 2 


v w «r « ^ I a. z a I J Z < V? w> ? 

I > i: Cj u r j w .? I i-u U £22 

l « 2 i. a O i j iii o 2 u “■* < « < 

s o < ojS) u! u y •- a; £ h <^>>> 


16 



Table 2 


Standard Deviations Assigned to the Observations 


Observation 

Standard Deviation 

GRARR R (range) 

10 meters 

R (range- rate) 

3 cm/sec 

Laser R 

1 meter 

Optical 6 (declination) 

2 seconds of arc 

(cos 8 )a (right ascension) 

2 seconds of arc 

Minitrack direction cosines 

3 xlflr 4 

NWL Doppler R 

4cm/sec 

Gravity Data, 

33/(N + 1)' A mgal 

5° x 5° mean gravity anomalies 


(N= no. of points given in 5° x 5° section) 
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GSFC Geopotential Solutions (Normalized Coefficients X 10°) 1 
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GSFC Geopotential Solutions (Normalized Coefficients X 10 6 ) 2 
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Table 4 


Station Coordinates 


LATITUDE 



3TN 

DDMM 

ss.ss 

ANCHOR 

14 

6117 

0.51 

LACRES 

103 

3216 

44.68 

LASHM2 

106 

5111 

9.30 

APLMND 

111 

39 9 

48.67 

1BP01N 

1021 

3825 

49.77 

1FTMYR 

1022 

2632 

53.45 

100MFR 

1024-3123 

25.19 

1SATAG 

1028-33 8 

58.49 

1MOJAV 

1030 

3519 

47.90 

1JOBUR 

1031-2553 

0 « 66 

1NEWFL 

1032 

4744 

29.84 

1GFORK 

1034 

48 1 

21.27 

1WNKFL 

1035 

5126 

46.48 

1ROSMN 

1037 

3512 

7.13 

lORORL 

1038-3537 

32.15 

1ROSMA 

1042 

3512 

7.33 

1 TAN AN 

1043-19 0 

31.72 

MADGAS 

1123-19 1 

14.19 

ROSRAN 

1126 

3511 

45.51 

ULASKR 

1128 

6458 

19.05 

PRETOR 

2115-2556 

48.21 

1UNDAK 

7034 

48 1 

21.41 

1EDINB 

7036 

2622 

46.76 

1COLBA 

7037 

3853 

36»l8 

1BERMD 

7039 

3221 

49.86 

1PURIO 

7040 

1815 

28.94 

1GSFCP 

7043 

39 1 

15.32 

1DENVR 

7045 

3938 

48.02 

GODLAS 

7050 

39 1 

14.53 

WALLAS 

7052 

3751 

36.19 

1JUM40 

7072 

27 1 

14.54 

1JUBC4 

7074 

27 l 

14.79 

1SUDBR 

7075 

4627 

20.40 

1 JAM AC 

7076 

18 4 

34.78 

1 ORGAN 

9001 

3225 

25.04 

10LFAN 

9002- 

2557 

35.95 

WOOMER 

9003-31 6 

1.98 

1SPAIN 

9004 

3627 

46.74 

1 TOKYO 

9005 

3540 

23.00 

1 NATAL 

9006 

2921 

34.71 

1QUIPA 

9007-1627 

56.74 

1SHRAZ 

9008 

2938 

13.87 

1CURAC 

9009 

12 5 

25.19 

1JUPTR 

9010 

27 1 

14.15 

1VILDO 

9011-3156 

34.67 

1MAUIO 

9012 

2042 

26.15 

AUSBAK 

9023-3123 

25.82 


LONGITUDE 

HEIGHT 

DDDMM 

SS.SS 

METERS 

21010 

29.69 

68.7 

25314 

45.42 

1149.0 

35858 

26.32 

221.5 

283 6 

12.78 

97.3 

28254 

48.71 

-40.1 

278 8 

4.33 

-31.1 

13652 

15.89 

129.9 

28919 

53.60 

707.2 

243 5 

59.30 

888.8 

2742 

26.59 

1539.0 

30716 

46.34 

69.0 

26259 

19.87 

219.4 

35918 

8.38 

104.0 

277 7 

41.75 

865.4 

14857 

14.80 

939.5 

27 7 7 

41.44 

865.5 

4717 

59.73 

1364.4 

4718 

11.75 

1390.8 

277 7 

26.20 

828.9 

21229 

12.44 

341.6 

282C 

51.88 

1579.8 

26259 

19.60 

217.3 

26140 

7.73 

21.4 

26 7*7 

41.11 

227*3 

2952C 

35.37 

-15.2 

294 C 

23.77 

-5.7 

2831C 

20.33 

0. 6 

25523 

38.87 

1759.3 

2831C 

18.33 

12.3 

28429 

23.76 

-42.1 

27953 

13.07 

-30.6 

27953 

13.11 

-34.2 

279 3 

10.99 

233.3 

28311 

27.12 

422.4 

25326 

49.06 

1616.4 

2814 

52.83 

1553.5 

13647 

3.63 

151.0 

35347 

37.13 

56.9 

13932 

16.65 

81.0 

7927 

27.60 

1869.6 

28830 

24.82 

2476.5 

5231 

11.52 

1578.2 

291 9 

44.71 

-28.8 

27953 

13.55 

-26.8 

29453 

36.93 

613.9 

20344 

33.98 

3041.1 

13652 

43.96 

128.4 
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Table 6 


Degree Variances for Gravity Anomalies (Mgal 2 ) 


DEGREE 

GSFC 

COMB. 

S.E.Z 

GSFC 

SAT. 

3 

33.577 

32.842 

33.344 

4 

21.509 

21.805 

21.611 

5 

20.359 

17.785 

19.790 

6 

19.651 

15.652 

18.764 

7 

17.649 

15.491 

19.741 

8 

9.029 

6.639 

10.311 

9 

10.538 

12.651 

11.454 

10 

9.951 

12.860 

11.158 

11 

7.374 

12.234 

8.174 

12 

4.415 

5.099 

4.877 

13 

13.375 

11.121 

5.001 

14 

11.333 

8.431 

2.424 

15 

11.204 

13*215 

2.301 

16 

6.096 

13.844 

0.630 

TOTAL 

196 . 

199 . 
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Table 7 

Comparison of Coefficients with S, E. II 

Normalized Coefficients Multiplied by 10 6 

RMS OF COEFFICIENT DIFFERENCES BY DEGREE 


DEGREE 

GSFC 

COMB. 

GSFC 

SAT. 

SEI 

B13.1 

3 

0*021 

0*032 

0*112 

0.068 

4 

0*039 

0*043 

0*064 

0.030 

5 

0*049 

0*059 

0*055 

0.056 

6 

0*062 

0*065 

0*072 

0.044 

7 

0 * 0 S 9 

0*073 

0*093 

0.060 

8 

0*031 

0*042 

0*065 

0.061 

9 

0*046 

0*050 

0*0 

0.051 

10 

0*052 

0*063 

o 

• 

o 

0.058 

11 

0 *066 

0*057 

0.0 

0.056 

12 

0*044 

0*048 

0.0 

0.050 

13 

0*065 

0*055 

0.0 

0.066 

14 

0*049 

0*023 

0.0 

0.060 

IS 

0.047 

0*025 

0.0 

0.048 

16 

0*034 

0*018 

0*0 

0.017 

TOTAL RMS 

0*050 

0*052 

0.077 

0.053 
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Table 8 


Station Geoid Heights from the Combination Solution 


and MSL Heights from Survey 


DATA 

STATION 

STATION 

LAT* 

LON* 

MSL 

GEOID 

TYPE 

NAME 

NUMBER 

(DEG) 

(DEG) 

(M) 

(«) 

0 

ANCHOR 

14 

61.28 

210.17 

68. 

7. 

0 

L4CRES 

103 

32.27 

253.23 

1203. 

-24. 

D 

LASHM2 

IC6 

51. 18 

358.97 

190. 

48. 

D 

APLMND 

111 

39.15 

283,10 

14i. 

-38, 

D 

PRETOR 

2115 

-25.93 

28.33 

1580. 

24. 

M 

18P0IN 

1021 

38.42 

282. 9C 

6. 

-38. 

M 

1FTMYR 

1022 

26.53 

278.13 

5. 

-27. 

M 

100MFR 

1024 

-31.38 

136,87 

133. 

-1. 

M 

1SATAG 

1028 

-33.13 

289.32 

693. 

26. 

M 

1M0JAV 

1030 

35.32 

243.08 

929. 

-30. 

M 

1JOBUR 

1031 

-25.88 

27.70 

1522. 

25. 


1NEWFL 

1032 

47.73 

307.27 

69. 

8. 


1GF0RK 

1034 

48.02 

262.98 

253. 

-26. 

M 

1WNKFL 

1035 

51.43 

359.30 

67. 

48. 

M 

1ROSMN 

103? 

35.20 

277.12 

909. 

-35. 

M 

10RORL 

1038 

-35.62 

148.95 

932. 

16. 

M 

1ROSMA 

1042 

35.20 

277.12 

904. 

-35. 


IT AN AN 

1043 

-19.00 

47.28 

1378. 

-7. 

M 

1UNDAK 

7034 

48.02 

262.98 

253. 

-26. 

V 

1E0IN8 

7036 

26.37 

261.67 

60. 

-19. 

P 

1COLRA 

7037 

38.88 

267.78 

273. 

-31. 

M 

1BERMD 

7039 

32.35 

295.33 

31. 

-45. 

M 

1PURI0 

7040 

18.25 

294. OC 

50. 

-46. 

M 

1GSFCP 

7043 

39.02 

283.17 

53. 

-38. 

M 

1CENVR 

7045 

39.63 

255.38 

1790. 

-20. 

M 

1JUM40 

7072 

27.02 

279.88 

14. 

-31. 

« 

ISUDBR 

7075 

46.45 

279.05 

281. 

-36. 

M 

1JAMAC 

7076 

18.07 

283.18 

446. 

-21. 

G 

MADGAS 

1123 

-19.02 

47. 3C 

1399. 

-7. 

G 

RGSRAN 

1126 

35.18 

277.12 

874. 

-33. 

G 

ULASKR 

1128 

64.97 

212.48 

349. 

rr 

1 • 

t 

GODLAS 

7050 

39.02 

283.17 

55. 

-38. 

L 

WALLAS 

7052 

37.85 

284.48 

9. 

-39. 

8 

10RGAN 

9001 

32.42 

253.43 

1651. 

-24. 

8 

10LFAN 

9002 

-25.95 

28.23 

1544. 

24. 

B 

1SPAIN 

9004 

36.45 

353.78 

26. 

51. 

B 

1TUKYU 

9005 

35.67 

139.53 

60. 

39. 

0 

1NATAL 

9006 

29.35 

79.45 

1927. 

-49. 

8 

JLQUIPA 

9007 

-16.45 

288.50 

2452. 

31. 

B 

ISHRAZ 

9008 

29.63 

52.52 

1596. 

-9. 

8 

ICURAC 

9009 

12.08 

291.15 

7. 

-25. 

B 

1JUPTR 

9010 

27.02 

279.08 

15. 

-31. 

B 

1VILOO 

9011 

-31.93 

294.88 

598. 

27. 

B 

1MAUIO 

9012 

20.70 

203.73 

3034, 

6. 

B 

AUSBAK 

9023 

-31.38 

136.87 

141. 

-2. 

B 

WOOMER 

9005 

-51 *58 

156.87 

162, 

-2. 


*See Table 4 for more complete values. 
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Table 9 

Long** Arc Tests on Mean Elements for Resonant Satellites 



INTELSAT* 

2-F1 

Cosmos 41* 
Rocket 

Cosmos 382** 
Rocket 

Beat Period (days) 

900 

900 

125 

Data Span (days) 

1100 

2000 

430 

Gravity Fields 

RMS 

Mean Anomaly 
(degs.) 

RMS 

Mean Anomaly 
(degs.) 

RMS 

Semi-Major Axis 
(meters) 

1. S.E. n 

2.0 

6.5 

55.5 

2. GSFC-Comb. 

0.5 

3.6 

15.6 


Resonant Coefficients 

*<2,2), (3,2), (4,2) (4,4), (5,4), (6,4) (6,6), (7,6) .... etc. 

**(9,9), (10,9), (11,9) .... etc. 
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Table 10 

Summary of Solution, Data, and Results 


A. SOLUTION 

Geopotential solution in spherical harmonic series (degree n, order m). 


Type Solution Complete 
(n x m) 


Zonals 

Complete 


Resonant Coefficients 


Satellite (GEM 1) 12x12 

Satellite & 

Gravimetry 

(GEM 2) 16 x 16 


21 m = 12, 13, 14, n = m to 22 

m = 9, n = 9 to 16 
m = 11, n = 11, 12, 14 

21 Complete as above 


46 Center of mass tracking station locations (see tracking systems below) 


B. DATA 

Satellite Data 


Tracking Systems 25 Satellites (Data Processed) 

13 Baker-Nunn Cameras 300 weekly arcs of optical data 

23 MOTS Cameras 15 weekly arcs of electronic/laser data 

3 GRARR (range/ range- rate) 66 one/two-day arcs of GEOS-I and II 
2 LASERS flashing light data 

5 DOPPLER 20 weekly arcs of Minitrack data 

Gravimetry Data 

19,000 1° x 1° mean gravity anomalies of ACIC data, and 
2,000 1° x 1° mean gravity anomalies from other sources, used to form 
1,705 5° x 5° means, providing 70% world coverage 

C. RESULTS 

1. Station heights indicate a e = 6378145 meters. 

2. Mean value of gravity (Ag 0 = 3.28 mgal) indicates a e = 6378142 (m). 

3. Solution for 327 zonal and tesseral coefficients including Ag 0 
from satellite and gravimetry data. 

4. Solution using optical and electronic tracking data for location 
of 46 tracking stations. 
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